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A countrywide survey in Oman revealed Middle East 
respiratory syndrome coronavirus (MERS-CoV) nucleic 
acid in five of 76 dromedary camels. Camel-derived 
MERS-CoV sequences (3,754 nucleotides assem- 
bled from partial sequences of the open reading 
frame (ORF)1ia, spike, and ORF4b genes) from Oman 
and Qatar were slightly different from each other, 
but closely related to human MERS-CoV sequences 
from the same geographical areas, suggesting local 
zoonotic transmission. High viral loads in nasal and 
conjunctival swabs suggest possible transmission by 
the respiratory route. 


Background 

In June 2012 a novel betacoronavirus, subsequently 
named Middle East respiratory syndrome coronavirus 
(MERS-CoV), was isolated from a patient with fever and 
respiratory symptoms who had been admitted to a hos- 
pital in Jeddah, Saudi Arabia [1]. As of 22 April 2014, 
the number of reported laboratory-confirmed cases of 
MERS-CoV worldwide amounts to 333, with 107 deaths 
[2]. To date, only fifteen cases, including ten in Europe, 
have been reported outside of the Middle East [2] and 
the vast majority of cases were reported from Saudi 
Arabia. Thus, it seems that the virus may originate 
from the Arabian Peninsula. Family, healthcare associ- 
ated and community case clusters of MERS-CoV infec- 
tions have been reported (e.g. [3]). 


Besides limited human-to-human transmission, how- 
ever, epidemiological data point towards an animal 
reservoir of MERS-CovV. First evidence of such a reser- 
voir host was provided when all 50 investigated sera of 
dromedary camels (Camelus dromedaries) from Oman 
exhibited high-titre neutralising antibodies against 
MERS-CoV [4]. This observation has meanwhile been 
confirmed by several studies from other countries on 
the Arabian Peninsula and beyond (e.g. Egypt [5]). Other 
important farm animals on the Arabian Peninsula such 
as cattle, goats, sheep and chickens were also inves- 
tigated but were found negative [6,7]. Although differ- 
ent species of bats carry a variety of coronaviruses and 
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have been suggested as the most likely primary animal 
reservoir for MERS-CoV, so far only a short (190 nucleo- 
tides (nt)) sequence in a conserved region of the MERS- 
CoV genome was amplified from a faecal pellet of a bat 
(Taphozous perforates) in Saudi Arabia [8]. The present 
study was initiated in order to identify the virus in cam- 
els and to compare it genetically with human-derived 
MERS-CoV. 


Sample collection 

In December 2013, nasal and conjunctival swabs 
were taken from 76 dromedary camels of different 
age, breed and sex from all over Oman. All swabs 
were taken in duplicates, one stored in virus isolation 
medium and the other in the virus-inactivating buffer 
DNA/RNA Shield (Zymo Research, Irvine, USA). The 
former samples were frozen at -80°C at the Veterinary 
Research Center of the Ministry of Agriculture and 
Fisheries, Oman, for later use, while the latter samples 
were shipped to the University of Veterinary Medicine 
Vienna for analysis. DNA/RNA Shield effectively lyses 
cells and inactivates nucleases and infectious agents. 
In addition it ensures nucleic acid stability during sam- 
ple storage and transport at ambient temperatures. 


Nucleic acid extraction and polymerase 
chain reactions 

RNA was extracted employing a Quick-RNA MiniPrep 
Kit (Zymo Research, Irvine, USA), following the manu- 
facturer’s instructions. 


For screening, two published MERS-CoV reverse 
transcription-quantitative polymerase chain reaction 
(RT-qPCR) assays were performed, one in the open 
reading frame (ORF)1a gene region [9] and the other in 
the ORF1b gene region [10], and both were optimised 
for SuperScript III Platinium One-Step RT-qPCR System 
(Invitrogen/LifeTechnologies, CA, USA). For confirma- 
tion, the positive samples were retested by reverse 
transcription-polymerase chain reactions (RT-PCRs) in 
the RNA-dependent RNA polymerase (RdRp)- and nucle- 
ocapsid (N) gene regions [9] as well as in the spike 


TABLE 


Primers used for amplification of Middle East respiratory syndrome coronavirus genetic sequences and subsequent 


phylogenetic analysis, 2013 


Primer name, position? and sequence 
MERS_1767_F: 5'-CTCGCAATTCTCTCTGGAAC-3' 


Target sequence Expected size 


MERS_2615_R: 5'-GTCAGTAGGTTGGAGCAGTC-3' Okita paint 
MERS_11419_F: 5'-CAA GCC CCA TTG CCT ATC 1G-3' 

MERS_12064_R: 5'-GCT TGA AGT ACG CTA GGA GTG-3' ORFia 645 nt 
MERS_22074_F: 5'-CGTAATGCCAGTCTGAACTC-3' eae ra 
MERS_23127_R: 5'-CAGGGTGAGTATTGATTAGCG-3’ P 95 
MERS_24156_F: 5'-GCTGATCCTGGTTATATGCAAGG-3' uke x 
MERS_24903_R: 5'-CAACCTCAATGTGGTTGCTAGG-3' P 747 
MERS_26042_F: 5'-CTT TGG CCA AAC AGG ACG CA-3' Saris aan 
MERS_26856_R: 5'-GAC GCC GAG AAA GCC ATA GTT C-3' 4 ‘i 


F: forward; MERS: Middle East respiratory syndrome; nt: nucleotide; ORF: open reading frame; R: reverse. 


@ Referring to the sequence of the MERS-coronavirus strain HCoV-EMC/2012 with the GenBank accession number: JX869059. 


gene region [11] and by five other RT-PCRs designed for 
subsequent phylogenetic analysis (Table). The MERS- 
CoV genome organisation and the location of the five 
RT-PCR amplicons used to generate a concatenated 
sequence for phylogenetic analysis are displayed in 
Figure 1. 


Allconventional RT-PCRs were conducted using OneStep 
RT-PCR kit (Qiagen, Hilden, Germany). Primer synthe- 
sis and sequencing in both directions were carried out 
by Microsynth (Balgach, Switzerland). The obtained 
sequences were verified by Basic Local Alignment 
Search Tool (BLAST) search, aligned using the ALIGN 
PLUS programme (Scientific & Educational Software), 
and compiled to one concatenated sequence. 


Phylogenetic analysis 

For phylogenetic analysis, three camel-derived (from 
Oman, Qatar, and Egypt) and 33 human-derived MERS- 
CoV sequences were included. Amultiple sequence align- 
ment was performed using BioEdit Sequence Alignment 
Editor version 7.0.9.0 and verified using the CLUSTAL_X 
programme (version 1.8). Phylogenetic neighbour-join- 
ing and maximum likelihood analyses were conducted 


FIGURE 1 


with the help of the Molecular Evolutionary Genetics 
Analysis (MEGA) 5 programme [12]. The evolutionary 
distances were computed using the Kimura 2-param- 
eter model. Bootstrap resampling analysis with 1,000 
replicates was employed. The sequences used for phy- 
logenetic analysis were deposited at GenBank under 
the accession numbers KJ573789-KJ573793 (MERS- 
CoV sequences derived from camel Oman_30_2013), 
and KJ598493-KJ598496 and KF933384 (sequences 
derived from camel Qatar_1_2013). 


Results 

Nasal and conjunctival swabs of five of a total of 76 
camels (6.6%) proved positive in all applied RT-qPCR 
and RT-PCR assays. The cycle threshold (Ct) values 
ranged from 15.74 to 36.29. Concatenated sequences 
of a total of 3,754 nts were obtained from all five posi- 
tive Omani camels by compiling the sequences derived 
from the five different PCR amplification products cov- 
ering the ORF1a, spike and ORF4b gene regions (Table, 
Figure 1). The five concatenated sequences derived 
from the Omani camels were 100% identical to each 
other, and they exhibited 99% identity (differing in 12 
of 3,754 nts) to the sequence derived from the Qatari 


Genome organisation of Middle East respiratory syndrome coronavirus and location of five reverse transcription- 
polymerase chain reaction amplification products used to generate a concatenated sequence for phylogenetic analysis, 2013 
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FIGURE 2 


Phylogenetic analysis of three camel- and 33 human-derived Middle East respiratory syndrome coronavirus (MERS-CoV) 
nucleotide sequences, 2013 


KF186567, Human, Saudi Arabia, Al-Hasa_1_2013 
KF186566, Human, Saudi Arabia, Al-Hasa_2_2013 
KF186565, Human, Saudi Arabia, Al-Hasa_3_ 2013 
KF600627, Human, Saudi Arabia, Al-Hasa_12_2013 
KF600634, Human, Saudi Arabia, Al-Hasa_21_2013 
KF600645, Human, Saudi Arabia, Al-Hasa_15_2013 
KF600651, Human, Saudi Arabia, Al-Hasa_18_2013 
KJ156866, Human, Saudi Arabia, Al-Hasa_25_ 2013 
KJ156869, Human, Saudi Arabia, Riyadh_9_2013 
KF186564, Human, Saudi Arabia, Al-Hasa_4_2013 
79 KF600647, Human, Saudi Arabia, Al-Hasa_17_2013 
KF600632, Human, Saudi Arabia, Al-Hasa_19_2013 
KF600652, Human, Saudi Arabia, Riyadh_2_2012 
KF745068, Human, France, France/UAE_2013 
@ KJ573789-93, Camel MERS-CoV/Oman_30_2013 
KF600630, Human, Saudi Arabia, Buraidah_1_2013 
84 KJ156910, Human, Saudi Arabia, Hafr-Al-Batin_2_2013 Clade B 
KJ156874, Human, Saudi Arabia, Hafr-Al-Batin_6_2013 
KF600628, Human, Saudi Arabia, Hafr-Al-Batin_1_2013 


. KJ156934, Human, Saudi Arabia, Riyadh_14_2013 
86 KF961222, Human, Qatar_4_2013 
@ KJ598493-96 & KF933384, Camel MERS-CoV/Qatar_1_2013 
KF961221, Human, Qatar_3_2013 
94 KF600620, Human, Saudi Arabia, Bisha_1_2012 


97 ' KF600612, Human, Saudi Arabia, Riyadh_1_2012 
99 KC164505, Human, UK, England_1_2012 
KC667074, Human, UK, England/Qatar_2012 
KF600613, Human, Saudi Arabia, Riyadh_3_2013 
- KJ156952, Human, Saudi Arabia, Riyadh_4_2013 
KF192507, Human, UAE, strain Munich_2013 
KJ156949, Human, Saudi Arabia, Taif_1_2013 
KJ156944, Human, Saudi Arabia, Riyadh_5_2013 
KJ156881, Human, Saudi Arabia, Wadi-Ad-Dawasir_1_2013 
® KJ477102, Camel NRCE-HKU205, Egypt, 2013 
JX869059, Human, Saudi Arabia, Jeddah, strain HCoV_EMC, 2012 CladeA 
KC776174, Human, Jordan, strain Jordan_N3/2012 


0.0005 


UAE: United Arab Emirates; UK: United Kingdom. 


Each 3,754 nucleotide long sequence used to generate the tree was obtained from concatenating partial sequences of the open reading frame 
(ORF)1a, spike and ORF4b gene regions. Of note the different clustering of the camel-derived sequences originating from Oman (marked 


with a red diamond), Qatar (blue diamond) and Egypt (green diamond). The Qatari and Omani camel-derived MERS-CoV sequences cluster 
close to the human-derived sequences originating from the same areas. 
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camel and also 99% identity (31 nts difference) to the 
known respective sequence derived from an Egyptian 
camel (Figure 2). 


The phylogenetic analysis involving all three currently 
available camel-derived MERS-CoV sequences and all 
33 available human-derived sequences is shown in 
Figure 2. It clearly indicates that the camel-derived 
MERS-CoV sequences are clustering independently 
from each other, but together with the human-derived 
MERS-CoV sequences from the same geographical 
areas. This could be demonstrated unambiguously 
for the Qatari samples (exhibiting only 1 nt difference 
between the camel-derived and the human-derived 
MERS-CoV sequences), but also for the camel-derived 
Omani sequence, which is clustering close to a human 
MERS-CoV sequence from neighbouring United Arab 
Emirates (GenBank accession number: KF745068, 
with only 5 of 3,754 nts being different; no human 
Omani MERS-CoV sequences are available so far). 
For the Egyptian camel-derived MERS-CoV sequence 
[13] no human-derived sequence with as high relative 
sequence identity and/or from a respectively close 
geographical area has yet been detected. 


Discussion 

The results of our study and similar studies from Qatar 
[11] and Saudi Arabia [14-16] show a close genetic rela- 
tionship between camel-derived and human-derived 
MERS-CoV from the same geographical areas, suggest- 
ing local zoonotic transmission. A proof of cross-species 
transmission of MERS-CoV from dromedary camels to 
humans was reported in a publication on human infec- 
tion with MERS-CoV after exposure to infected camels 
in Saudi Arabia recently [15]. The authors concluded 
that camels may act as a direct source of human MERS- 
CoV infection. For the implementation of effective, but 
reasonable precautions, however, further pieces of the 
puzzle must be put together: a Ct-value of 15.74, cor- 
responding to approximately 33 million RNA copies, in 
the nasal swab of one of the Omani camels (this study) 
indicates a very high virus load and the probability of 
the transmission of high numbers of virus particles by 
the nasal route. Data from another recent publication 
suggest that neither passively acquired maternal anti- 
bodies nor prior infections seem to lead to full immu- 
nity in the camels, and viral nucleic acid was detected 
in such animals, though at a significantly lower level 
[16]; these and other recent results, however, require 
confirmation. 


Evidence was provided that MERS-CoV has been circu- 
lating in camels at least since 1992 [14], but probably 
much longer, while strains in humans emerged only 
recently [17]. Initially, dromedary camels have been 
considered possible intermediate hosts of MERS-CoV, 
while bats were suggested as primary reservoirs — in 
analogy to severe acute respiratory syndrome-corona- 
virus (SARS-CoV). According to our current knowledge, 
however, dromedary camels might even be considered 
the primary reservoir hosts of MERS-CoV. Comparative 


analyses of camel- and human-derived MERS-CoV 
sequences exhibit »99% nt identities to each other with 
so far no striking genetic differences observed. Given 
the high prevalence of MERS-CoV-specific antibodies in 
dromedary camels [4,5], infection of camels with this 
virus seems to be frequent on the Arabian Peninsula. 
This is in sharp contrast to the comparatively low num- 
ber of recorded human infections. Since currently no 
other obvious explanation for this discrepancy exists, 
it is postulated that a high infectious dose through very 
close contact between an infected camel and a human 
being is required for initiation of human MERS-CoV 
infection by camels. Younger camels may play a par- 
ticular role in zoonotic transmission since they seem 
to be more frequently infected and seem to shed more 
virus than older ones [14-16]. Although the respiratory 
route is, in our opinion, the most likely route of trans- 
mission, a recent paper demonstrated that MERS-CoV 
survived in raw camel milk slightly longer than in milk 
of other species [18], suggesting further investigations 
on a possible food-borne route of transmission. Also, it 
is not clear yet whether the course of MERS-CovV infec- 
tion in camels is generally asymptomatic or associated 
with mild respiratory symptoms, as suggested recently 
[16]. 


Conclusions 

Phylogenetic analysis and high MERS-CovV viral load 
in nasal swabs of dromedary camels suggest local 
zoonotic transmission through the respiratory route. 
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